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On the basis of perturbative QCD and the relativistic quark model we calculate 
relativistic and bound state corrections in the production processes of a pair of S-wave 
and P-wave charmonium states. Relativistic factors in the production amplitude 
connected with the relative motion of heavy quarks and the transformation law 
of the bound state wave function to the reference frame of the moving S- and P- 
wave mesons are taken into account. For the gluon and quark propagators entering 
the production vertex function we use a truncated expansion in the ratio of the 
relative quark momenta to the center-of-mass energy ^/s up to the second order. 
The exact relativistic treatment of the wave functions makes all such second order 
terms convergent, thus allowing the reliable calculation of their contributions to the 
production cross section. Relativistic corrections to the quark bound state wave 
functions in the rest frame are considered by means of the Breit-like potential. It 
turns out that the examined effects change essentially the nonrelativistic results of 
the cross section for the reaction e"*" + e~ — > J/^'(r/c) + Xcj(^c) at the center-of-mass 
energy ^/s = 10.6 GeV. 
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I. INTRODUCTION 

The large rate for the exclusive double charmonium production measured at the Belle 
and BaBar experiments [H, 0] reveals definite problems in the theoretical description of these 
processes [sl 0, • Many theoretical efforts were made in order to improve the calculation 
of the production cross section + e" J 1^ + rjc- They included the analysis of other 
production mechanisms for the state J I'^+rjc p,0] and the calculation of different corrections 
which could change essentially the initial nonrelativistic result js], [o], [l^, |Tl|, 12, 13, 14, 15, 
16, [O, Despite the evident successes achieved on the basis of NRQCD, the light cone 
method, quark potential models in order to resolve the discrepancy between the theory and 
experiment, the double charmonium production in e^e~ annihilation remains an interesting 
task. On the one hand, the reason is that there exist the production processes of the P- 
and D-wave charmonium states which should be investigated as the production of S-wave 
states. On the other hand, the variety of the used approaches and the model parameters in 
this problem raises the question about the comparison of the obtained results resulting in 
a better understanding of the quark-gluon dynamics. Two sources of the enhancement of 
the nonrelativistic cross section for the double charmonium production are revealed to the 
present: the radiative corrections of order 0{as) and relative motion of c-quarks forming 
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FIG. 1: The production amplitude of a pair of S- and P-wave charmonium states in e~^e~ annihi- 
lation. S denotes the S-state meson and V the P-wave meson. The wavy line shows the virtual 
photon and the dashed line corresponds to the gluon. T is the production vertex function. 



the bound states. 

In this work we continue the investigation of the exclusive double charmonium production 
in e~^e~ annihilation on the basis of a relativistic quark model 

HI 

Il9l . |20j, l2l|] in the case of S- 

and P-wave charmonium states. The relativistic quark model provides the solution in many 
tasks of heavy quark physics. In particular, it gives the possibility to study the question 
about a broadening of the meson wave functions due to the account of special corrections 
which can lead to the increase of the double charmonium production cross sections. Thus, 
the aim of this study consists in the calculation of the relativistic effects in the processes 
e"*" + e~ J/^ij]^ + yrjjhr) on the basis of a relativistic approach to the quarkonium 
production suggested in Refs. [13, 1^ • 



II. GENERAL FORMALISM 



We consider the following reactions e"*" + J j'^ijjc) + Xcj{hc), where the final state 

consists of the pair of S-wave (J/^ or rjc) and P-wave {Xco, Xa, Xc2 or he) charm mesons. 
The diagrams that give contributions to the amplitude of these processes in the leading 
order of the QCD coupling constant as are presented in Fig.l. Two other diagrams can be 
obtained by corresponding permutations. There are two stages of the production process. 
In the first stage, which is described by perturbative QCD, the virtual photon 7* produces 
four heavy c-quarks and c-antiquarks with the following four-momenta: 

Pi,2 = ^P±P, {p-P) = 0; q,,2 = lQ±q, {q ■ Q) = 0, (1) 

where P{Q) are the total four-momenta, p = Lp(0,p), q = Lp(0, q) are the relative four- 
momenta obtained from the rest frame four-momenta (0, p) and (0, q) by the Lorentz trans- 
formation to the system moving with the momenta P,Q. In the second nonperturbative 
stage, quark-antiquark pairs form the final mesons. 

Let consider the production amplitude of the S-wave vector state [J/"^) and P-wave 
states Xcj (J=0,l,2), which can be presented in the form [3, [2l|, 22, [2^: 



xSp {^^{p, P)rf (p, q, P, g)vl/^(g, g)7. + ^""{q, g)rf (p, q, P, Q)^^ip, P)7.} , 

where is the QCD coupling constant, a is the fine structure constant and Qc is the c 
quark electric charge. The relativistic S- and P-wave functions of the bound quarks '^'^'^ 
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accounting for the transformation from the rest frame to the moving one with four momenta 
P, Q, are 
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where Vi = P/Mg, V2 = Q/M-p; es is the polarization vector of the vector charmonium J/\I/; 
ep{Q,Sz) is the polarization vector of the spin-triplet state Xcj, ^{p) = y/p^'-Tm? and m 
is the c quark mass. At the leading order in as the vertex functions 2(p, P; Q, Q) can be 
written as 



(/ — gi)^ — + ie 
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where the gluon momenta are ki = pi + qi, k2 = P2 + ^2 and = s = (P + Q)^ = (p_ +P+Y1 
P-, p+ are four momenta of the electron and positron. The dependence on the relative 
momenta of c-quarks is presented both in the gluon propagator D^y{k) and quark propagator 
as well as in the relativistic wave functions (3), (4). Taking into account that the ratio of 
the relative quark momenta p and q to the energy ^/s is small, we expand the inverse 
denominators of quark and gluon propagators as follows: 
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(7) 



(8) 



(9) 



In the expansions (7)- (9) we accounted for terms of second order in relative momenta p 
and third order in relative momenta q. Substituting (7)-(9), (3)-(4) in (2) we preserve 
relativistic factors entering the denominators of the relativistic wave functions (3)- (4), but 
in the numerator of the amplitude (2) we take into account corrections of second order 
in |p|/m and up to fourth order in |q|/m. This provides the convergence of the resulting 
momentum integrals. Then the angular integrals are calculated using the following relations: 



{g,jiv - vi ^vi ^) p'^Rs{p)dp, 



(10) 
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Qc.qgqj'^oifijj^ = -^;^[£^{Q, L,)Pc,p + e^iQ, L,)P^p + ep{Q, L,)Pa^] q^Rp{q)dq, 

(12) 

where Pap = {gap — V2 0^2 /?), Rs{p), Rv{q) are the radial momentum wave functions of 
S- and P-wave charmonium states, e^{Q,Lz) is the polarization vector in orbital space. 
The integrals in (10) and (12) look formally divergent, but the original momentum integrals 
contain also definite relativistic factors which lead to their convergence. For a specific P-wave 
state, summing over and in the amplitude (2) can be further simplified as ^3] 



-^{9ap-V2 aV2 J = 0, 

^ <1,L,;1,S,|J,J, JQ,L,)£;^(g,5,) = <' -i=e,^,,t;2-e*^(g,J,), J = 1, (13) 
^""-^ Ua/3(Q,^.), ^ = 2, 

where < 1, Lz]l, Sz\J, Jz > are the Clebsch-Gordon coefficients. Calculating the trace in 
the amplitude (2) by means of expressions (3)-(6), (13) and the system FORM 25|], we find 
that the tensor parts of the four amplitudes describing the production of S- and P-wave 
charmonium states in the used approximation have the following structure: 

Si,f3{J/^ + Xco) = Aie*s p + A2V1 p{v2e*s) + ^3^72 /3(^24). (14) 

S2AJI'^ + Xci) = B,eaXyf3V?V2^* HQ, Jz){v2e*s) + ^s^^a^/S^^ HQ^ Jz) + (15) 
+BsVi /BEax-yaViV^Es^e* "{Q, Jz) + BiV2 uSax^^v'^v^ey's* ''{Q, Jz), 
^3,/3(J/^+Xc2) = <.(g, J.)[Cie*s''g.f3+C2vnv2e*s)9p.+C3v'^e*s''v, ^+C4<£r^2 ,3+ (16) 
+C,vy^,e*s p + C6«^f (24) + Cjvyiv2 p{v2e*s)l 
S^AVc + K) = Divi p{vie*{Q, Lz)) + D2V2 p{vie*{Q, Lz)) + D^e*p{Q, Lz), (17) 

where the coefficients Ai, Bi, Ci, Di can be presented as sums of terms containing the 
factors u = Mp/{Mp + Ms), n = m/{Mp + Ms) and dj = c\p)c^{q) = [{m - e{p))/{m + 
e{p))Y[{m — e{q))/{m + e{q))]' , preserving terms with i+ j < 2, and = (Mp -|- Ms)"^ / s up 
to terms of order O(r^). Exact analytical expressions for these coefficients are sufficiently 
lengthy (compare with the results written in Appendix A of our previous paper [3]), so, 
we present in Appendix A of this work only their approximate numerical form using the 
observed meson masses and the c-quark mass m = 1.55 GeV. 

Introducing the scattering angle 6 between the electron momentum pe and the momentum 
P of the J/\l/ meson, we can calculate the differential cross section da /d cos 6 and then the 
total cross section cr as a function of r^. We find it useful to present the charmonium 
production cross sections in the following form {k = 0,1,2,3 corresponds to Xco, Xa, Xc2 
and he): 



d912k"='U^(1 — u)^ s[M'p + Ms) ~q 

(18) 
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where the functions Pj^^^ (k=0, 1,2,3) are written exphcitly in Appendix B, 
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The parameters Ui can be expressed in terms of momentum integrals /„, J„ as follows: 
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On the one side, in the potential quark model the relativistic corrections, connected 
with the relative motion of heavy c-quarks, enter the production amplitude (2) and the 
cross section (18) through the different relativistic factors. They are determined in the final 
expression (18) by the specific parameters cUj. The momentum integrals which determine the 
parameters Ui are convergent and we calculate them numerically, using the wave functions 
obtained by the numerical solution of the Schrodinger equation. The exact form of the 
wave functions ^E'q (p) and ^E'o'(q) is important for improving the accuracy of the calculation 
of the relativistic effects. It is sufficient to note that the double charmonium production 
cross section a{s) in the nonrelativistic approximation contains the factor |i?5(0)p|i?p(0)p. 
Small changes of the numerical values of the bound state wave functions at the origin lead 
to substantial changes of the final results. In the approach based on nonrelativistic QCD 
this problem is closely related to the determination of the color-singlet matrix elements for 
the charmonium j26|. Thus, on the other side, there are relativistic corrections to the bound 
state wave functions \l/g (p), \E'^(q). In order to take them into account, we suppose that 
the dynamics of a cc-pair is determined by the QCD generalization of the standard Breit 
Hamiltonian S,^,!^: 



H = Ho + AUi + Af/2, Ho = 2Vp2 + - 2m - 



+ Ar + B, 



(23) 
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where n/ is the number of flavors, Ca = 3 and Cp = 4/3 are the color factors of the SU(3) 
color group. For the dependence of the QCD coupling constant as(/x^) on the renormalization 
point /i^ we use the leading order result 



/3oln(/iVA')' 



(26) 
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TABLE I: Numerical values of the relativistic parameters (19), (20), (22) in the double charmonium 
production cross section (18). 



Meson (cc) 






^5(0), Gey3/2 




LOi{S) or (jJ4{V) 


U!2{S) or U!5{V) 








0.81 




-0.20 


0.0078 


Vc 


I' So 


0-+ 


0.92 




-0.20 


0.0087 


XcO 


l3Po 


0++ 




0.19 


-0.15 


0.0065 


Xcl 


I^PI 


1++ 




0.18 


-0.14 


0.0065 


Xc2 




2++ 




0.18 


-0.15 


0.0065 


he 




1+- 




0.18 


-0.14 


0.0065 



The typical momentum transfer scale in a quarkonium is of order of the quark mass, so we 
set the renormalization scale /i = m and A = 0.168 GeV, which gives as = 0.314 for the 
charmonium states. The parameters of the linear potential A = 0.18 GeV^ and B = —0.16 
GeV have usual values of quark models. Starting with the Hamiltonian (23) we construct the 
effective potential model based on the Schrodinger equation and find its numerical solutions 
in the case of S- and P-wave charmonium |3fl] . The details of the used model are presented 
in Appendix C. Then we calculate the matrix elements entering in the expressions for the 
parameters Ui and obtain the value of the production cross sections at •\/i=10.6 GeV. Basic 
parameters which determine our numerical results are collected in Table I. The comparison 
of the obtained results with the previous calculations 0, 0, [U, 11] and experimental data 
(ll, [2I is presented in Table II. 



III. NUMERICAL RESULTS AND DISCUSSION 

In this paper we have investigated the role of relativistic effects in the production processes 
of S- and P-wave mesons (cc) in the quark model. In the present study of the production 
amplitude (2) we kept relativistic corrections of two types. The first type is determined by 
several functions depending on the relative quark momenta p and q arising from the gluon 
propagator, the quark propagator and the relativistic meson wave functions. The second type 
of corrections originates from the perturbative treatment of the quark-antiquark interaction 
operator which leads to the different wave functions ^l/o(p) ^ind ^I/o'(q) for the S-wave and 
P-wave charmonium states, respectively. In addition, we systematically accounted for the 
bound state corrections working with the observed masses of S-wave mesons (J/^I/, f]c) and 
P-wave mesons {Xcj, he)- The calculated masses of S-wave and P-wave charmonium states 



agree well with experimental values [33[] (see Table III). Note that the basic parameters of 



the model are kept fixed from the previous calculations of the meson mass spectra and decay 



widths [2l|, |3J, |35j . The strong coupling constant entering the production amplitude (2) is 
taken to be 0^=0.24 in accordance with the relation (26) at = 2m. 

Numerical results and their comparison with several previous calculations and experimen- 
tal data are presented in Table II. Theoretically, there were two studies of the production 
J/^{r]c) + Xcj{K) in e+e" annihilation in NRQCD 0, i]. They give 2.4 fb and 6.7 fb for 
the production of J /'^ + Xco- Such spread in results is explained by the different numerical 
values of the used parameters, i.e. the matrix elements, the mass of c-quark m and the 
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strong coupling constant a^. The third investigation of the process J/\l/ + Xco was done in 
the hght-front formahsm in 3l| where the result 14.4 fb was obtained. The essential growth 



of the cross section in jSli at ^/s = 10.6 GeV is connected with the use of specific light cone 
wave functions describing the relative motion of heavy c-quarks. The fourth study of the 
reaction + e~ J/\E' + Xco was devoted to the next-to-leading order QCD corrections 
Here it was shown that a sharp increase of the production cross section (cr = 17.9 fb) 
can be derived with the account of NLO in contributions. 

The exclusive double charmonium production cross section presented in the form (18) 
is convenient for a comparison with the results of NRQCD. Indeed, in the nonrelativistic 
limit, when u = 1/2, k = 1/4, uJi = {i > 1)), = 16m^/s, the cross section (18) coincides 
with the calculation in [3|. In this limit the functions FQ^''(r^) transform into corresponding 
functions Fk from [3|. When we take into account bound state corrections working with 
observed meson masses, we get u = Mp/{Mp + Mg) ^ 1/2, h, = m/{Mp + Ms) ^ 1/4. This 
leads to the modification of the general factor in (18) in comparison with the nonrelativistic 
theory and the form of the functions Fq^-* (see |3|). It follows from the numerical values of 
the parameters cUj, presented in Table I, that the relativistic corrections amount to 15-^20% 
in the production amplitude. Moreover, the relativistic effects decrease the values of the 
parameters Rs{0), -Rp(O), which transform into Rs{0), -Rp(O). In all considered reactions 

+ e~ — >• J/'if{ric) + Xcj{hc) the relativistic effects increase the nonrelativistic cross section, 
but in the case of the production rjc + he the sum of bound state plus relativistic corrections 
decreases the nonrelativistic cross section. It is necessary to point out once again that the 
essential effect on the value of the production cross sections J/\l/(?7c) + Xcj{hc) belongs to 
the parameters RsiS^), -Rp(O)' Small changes in their values can lead to significant 

changes for the production cross sections. Comparing the values of the parameters Rs{0) 
(|i?j/*,„,(0)P = 0.9; 1.2 GeV^/^), R!p{0) (|i?p(0)|2 = 0.043 GeV^/^) obtained in this study on 
the basis of quark model and in [32)] we see that the values of the radial wave functions at the 
origin are very close, but our value for the derivative of the radial wave function at the origin 
is slightly smaller. The calculation of radiative corrections Ofa^) to the nonrelativistic cross 
section of the production J/^ + Xco was done recently in [3^. It evidently shows that one- 
loop corrections are considerable (factor K = 2.8 to nonrelativistic result). As a result, the 
total value of the cross section J/\l/ + Xco significantly increases. It should be noted that the 
difference between the theory and both BaBar and Belle experiments became threatening. 

We presented a systematic treatment of relativistic effects in the S- and P-wave double 
charmonium production in e~^e~ annihilation. We explicitly separated two different types of 
relativistic contributions to the production amplitudes. The first type includes the relativis- 
tic v/c corrections to the wave functions and their relativistic transformations. The second 
type includes the relativistic p/ ^/s corrections emerging from the expansion of the quark 
and gluon propagators. The latter corrections were taken into account up to the second 
order. It is important to note that the expansion parameter p/ -y/s is very small. In our anal- 
ysis of the production amplitudes we correctly take into account relativistic contributions 
of order 0(f^/c^) for the S-wave meson and corrections of orders 0(f^/c^) and 0(f^/c^) for 
the P-wave mesons. We cannot keep corrections of order 0{v'^/c^) for the S-wave part of 
the amplitude (2) because they become divergent if we use expansions (7)-(9). Therefore 
the basic theoretical uncertainty of our calculation is connected with the omitted terms of 
order 0(p^/m^). Taking into account that the average value of the heavy quark velocity 
squared in the charmonium is < >= 0.3, we expect that they should not exceed 30% of 
the obtained relativistic contribution. These theoretical errors in the calculated production 
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TABLE II: Comparison of the obtained results with previous theoretical predictions and experi- 
mental data. 



State 
H1H2 


<^BaBar X 
BrH2^charged>2 


BrH2^charged>2 

ifb) fl] 


(^NRQCD 
ifb) [3] 


^(fb) 

m 


[3iJ 


[32] 


Our result 

m 


J/"^ + XcO 


10.3±2.5l|:^ 


6.4 ± 1.7 ± 1.0 


2.40 ± 1.02 


6.7 


14.4 


17.9(6.35) 


4.79 ± 0.80 


J/^ + Xcl 






0.38 ±0.12 


1.1 






1.07 ±0.23 


J/"^ + Xc2 






0.69 ±0.13 


1.6 






1.10 ±0.13 


Vc + he 






0.308 ±0.017 








0.24 ±0.02 



cross section at y/s = 10.6 GeV are shown directly in Table II. We have neglected the terms 
in the cross section (18) containing the product of In and Jn with summary index > 2 be- 
cause their contribution has been found negligibly small. There are no another comparable 
uncertainties related to the choice of m or any other parameters of the model, since their 
values were fixed from our previous consideration of meson and baryon properties j2l|, 
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APPENDIX A: THE COEFFICIENTS A„ B,, d, D, ENTERING IN THE 
PRODUCTION AMPLITUDES (14)-(17) 

These coefficients are the sums of the terms containing the parameters 
u = Mp/{Mp + Ms) and k = m/{Mj, + Ms). We present Ai, Bi, d, Di in numer- 
ical form using the observed meson masses and the mass of c-quark m = 1.55 GeV. 

e+ + e- ^ J/^ + XcO 

18 55 18 55 
Ai = -7.05 ± + O.OlSr^ ± C2o(7.05 — + O.OlSr^) + (Al) 



c; 41 09 
±Co2(5.42 - - O.OlSr^) + Cio(68.16 — - U.lbr^ ± 0.089r^)± 

13 71 41 48 

±Coi(23.29 8.11r2 + 0.066r^) ± Cii{-1Q.^?> + + 62.06r2 - 12.66r^), 

ry A ry ^ 

A2 = -8.35 ± 2Mr^ + C2o(8.35 - 2.34r2) ± ^02(1. 80 - 2Mr^) + (A2) 
±Cio(30.38-26.93r2 + 4.00r^)±Coi(5.96-9.62r2 + 2.69r^)±Cn(-19.22±36.18r2-23.80r^), 

A3 = -0.99 ± 0.99r2 ± C2o(0.99 - 0.99r^) + ^02(1. 98 - 0.99r^) + (A3) 
±Cio(1.65 - 9.72r^ ± 1.50/) ± Coi(0.97 - 2.21r^ ± 0.94/) ± Cii(-1.62 ± 4:.02r^ - 5.24/). 
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Bi = -0.002 + 2.66r2 + C2o(0.002 - 2mr^) - 2mr^Co2 + (A4) 
+Cio(-6.17-14.23r2+4.59r^)+Coi(-4.62-7.14r2 + 3.13r^)+Cii(11.56+28.13r2-27.18r^), 

B, = -2.10 - ^ + C2o(2.10 + ^) + Cio(14.99 - ^ - 3.34r2 + 0.016r^) + (A5) 

9 28 23 21 

+Coi(7.29 - - 0.83r^ + O.Olr^) + Cu{-2AAA + + 10.13r^ - 0.98/), 

S3 = -1.63r2 + 1.63r2C2o + 1.63r^Co2 + Cio(10.64r2 - 2.87r^) + (A6) 
+Coi(5.43r2 - 1.20r^) + Cii(-16.41r2 + 13.44r^), 
= -0.002 + 0.82r2 + C2o(0.002 - 0.82r2) - 0.82r2Co2 + (A7) 
+Cio(-6.17 - 1.23r^ + 1.77r^) + Cn(11.56 + 1.35r2 - 11.46r^). 

e+ + e- ^ J/^ + Xc2 

Ci = -1.84 + 1.84C20 + Cio(7.50 - 1.87r2 + O.Olr^) + (A8) 

+Coi(3.72 - 0.97r^ + 0.02r^) + Cii(-23.15 + 11.46r2 - 0.95r^), 

C2 = 2.63r^ - 2.63r^C2o - 2mr'^Co2 + Cio(-13.28r2 + 4.55r^) + (A9) 

+Coi(-7.72r^ + 2.96r^) + Cii(44.33r2 - 30.95r^), 

C3 = 1.59r^ - 1.59r2C2o - 1.59r^Co2 + Cio(-5.58r^ + 2.85r^) + (AlO) 

+Coi(-3.91r2 + l.lOr^) + Cii(12.99r=^ - 14.29r^), 

= -0.80r2 - O.SOr^Cso + 0.80r2Co2 + Cio(1.41r2 - 1.75r^) + (All) 

+Coi(4.02r2 - 1.51/) + Cii(-4.43/ + 12.19/), 

C5 = -2.39/ - 2.39/C20 + 2.39/C02 + Cio(6.9l/ - 4.33/) + (A12) 

+Coi(10.29/ - 2.98/) + Cii(-32.09/ + 32.60/), 

Ce = -2.37/Cio + Cii/(-12.59/ + 1.65/), C7 = 3.29/Cio + Cii/(-12.59/ + 1.65/). 

(A13) 

e+ + (:~ —f ijc + he 

Di = -1.58 + 1.29r2 + C2o(l-58 - 1.29/) + Co2(l-58 - 1.29/) + (A14) 
+Cio(0.70 - 9.33/ + 2.97/) + Coi(l-35 - 5.29/ + l.Ol/) + Cii(-0.60 + 19.42/ - 16.51/), 
D2 = 0.14 - 1.14/ + C2o(-0.14 + 1.14/) + Co2(-3.73 + 1.14/) + (A15) 
+Cio(7.11 + 2.77/ + 2.37/) + Coi(-0.12 + 4.72/ - 1.54/) + Cii(-6.08 - 3.69/ - 12.92/), 

Ds = -1.49 + ^ + C2o(1.49 - ^) + Cio(12.21 - ^ - 2.35/ + 0.02/) + (A16) 

2 97 13 GT 

+Coi(6.76 - -— - 0.80/ + 0.02/) + Cn(-32.42 + — — + 13.22/ - 1.35/). 
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APPENDIX B: THE FUNCTIONS F^ir^) (K=0,l,2,3) ENTERING IN THE 



PRODUCTION CROSS SECTION (18) 

e+ + e- ^ J/^ + Xco 

Fo^°^ = 12.94r^ + 690.48r^ - 391.86r^ + 39.82r^ + 6.36r^°, (Bl) 

f[^^ = -43.18r^ - 4554.04r^ + 5941. 64r^ - 1942.48r^ + 89.61r^°, (B2) 

) = -25.87r2 - 1380.96r^ + 783. 71r^ - 79.64r^ - 12.72r^°, (B3) 

^ = 36.04r2 + 7634.68r^ - 15478.8r^ + 10011.2r^ - 2095.47r^°, (B4) 

Ff^ = -25.29r^ - 996.50r^ + 1961.67r^ - 957.35r^ + 97.10r^°, (B5) 

= -38.81r^ - 447.77r^ + 601. llr^ - 101.74r^ - 12.87r^°, (B6) 

Fg^^' = 12.36r^ + 349.41r^ - 1216.3r^ + 1370.5r*^ - 579.05r^°, (B7) 

Fy^^^ = 84.44r^ + 6382.47r^ - 15011.8r^ + 13175.8r^ - 4056.79r^°. (B8) 

F^^^ = 165.06r'^ - 248.34r^ + 75.88r® + 17.33r^°, (B9) 

F^^^ = -1655.33/ + 3347.61r^ - 1829.98r^ + 80.46r^°, (BIO) 

fP = -330.12r^ + 496.69r^ - 151.76r^ - 34.67r^°, (Bll) 

F^^' = 4263.79r^ - 10572.2r^ + 8306.62r^ - 1953.72r^°, (B12) 

^ = -752.14/ + 1660.5/ - 979.66/ + 60.00r^°, (B13) 

F^ = -237.01/ + 390.21/ - 117.85/ - 35.64r^°, (B14) 

F^^^ = 920.55r/ - 2520.4/ + 2230.78/ - 612.50r^°, (B15) 

F^^^ = 7012.32/ - 18599.5/ + 15704.7/ - 4122.93r^°. (B16) 
e+ + e- ^ J/^ + Xc2 

Fo^^^ = 23.41/ + 99.63/ - 247.93/ + 110.30/ + 27.28r^°, (B17) 

Ff ^ = -78.15r^ - 931.82/ + 2499.9/ - 1804.84/ + 137.34r^°, (B18) 

fP = -46.82r^ - 199.27/ + 495.86/ - 220.60/ - 54.56r^°, (B19) 

F^ = 65.23/ + 1952.35/ - 5759.03/ + 5433.52/ - 1541.06r^°, (B20) 

fP = -94.75/ - 1014.54/ + 2475.82/ - 1509.41/ + 50.08r^°, (B21) 

fP = -70.23r2 - 213.25/ + 608.04/ - 269.46r^ - 55.66r^°, (B22) 

F^ = 95.89r2 + 2191. 11/ - 5632.06/ + 4350.3/ - 954.85rr^°, (B23) 

Fi^'^ = 316.34/ + 7972.42/ - 22933.8/ + 20774.7/ - 5658.04r^°. (B24) 
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e+ + e" 



<3) 



(3) 



= 11.69r2 - 27.29r^ + 35.00r^ - 22.48r^ + 6.02r^°, 
-lO.SSr^ + 41.40r^ - 185.00r^ + 250.24r® - 141.63r^°, 



.(3) 



-23.375r^ + 54.57r^ - 70.00r^ + 44.97r*^ - 12.04r 
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.(3) 
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.(3) 



.(3) 



2.29r2 + 88.53r^ + 5.76r^ - 348.17r* + 451.27r 



-19.99r^ + 118.20r^ - 215.81r^ + 199.84r** - 86.92r 



10 



10 



.(3) 



.(3) 



-35.06r^ + 101. 73r^ - 115.65r^ + 61.33r*^ - 12.52r 
= 8.55r2 - 81.12r^ + 268.40r^ - 363.38r^ + 234.48r^°, 
17.70r2 - 263.94r^ + 959.68r^ - 1671. 7r*^ + 1317.3ri°. 



(B25) 
(B26) 
(B27) 
(B28) 
(B29) 
(B30) 
(B31) 
(B32) 



APPENDIX C: EFFECTIVE RELATIVISTIC HAMILTONIAN 

For the calculation of the relativistic corrections in the bound state wave functions , 
^'o') we consider the Breit potential (23). It contains a number of terms which should be 
Tansformed in order to use the program of numerical solution of the Schrodinger equation 
30 1 . The rationalization of the kinetic energy operator can be done in the following form 

(CI) 



T 



2\/ + 



where rh is the effective mass of heavy quarks, 



+ rri^ p^ 2m 
2 ^ ^ + — 

m E 



E 

m = — 
2 



(C2) 



Pgjj should be considered as a new parameter which effectively accounts for relativistic 
corrections in (CI). Numerical values of vlff fo^' S- and P-wave charmonium states discussed 
13l . Il8l . |20[ | are presented in Table III. In the case of S-wave states it is necessary to 



m 



transform the (5-like terms of the potential. For this aim, we use the known smeared 5- 
function of the Gaussian form [37] : 



5{v) 



TC 



3/2 



(C3) 



with the additional parameter b which defines the hyperfine splitting in the (cc) system. 
Since the numerical results are practically not dependent on b in the range of commonly 
used values 1.5 -j- 2.2, we take b = 1.5 GeV. The second term in the Breit potential (23), 
which also has to be transformed, takes the form: 



AU 



2as 



(C4) 
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TABLE III: The parameters of the effective relativistic Hamiltonian. 









771 (tPV 

1 1 L ^ V_.I V . V 


E GeV 


P, GeV 


b, GeV 


M*'^, GeV 








0.5 


0.85 


0.087 


0.75 


1.5 


3.044 


3.097 


Vc 




0.5 


0.85 


0.087 


0.75 


1.5 


2.989 


2.980 


XcO 


l3Po 


0.6 


0.87 


0.479 


0.55 




3.437 


3.415 


Xcl 


l^Pl 


0.6 


0.87 


0.479 


0.55 




3.479 


3.511 


Xc2 


1^P2 


0.6 


0.87 


0.479 


0.55 




3.520 


3.556 


he 




0.6 


0.87 


0.479 


0.55 




3.486 


3.526 



In order to replace it by the effective term containing the power-like potentials, we use the 
approximate charmonium wave functions which can be written for S- and P-wave states as 



3/4 




The wave functions (C5) give a good approximation of the true quark bound state wave 
functions in the region of nonrelativistic momenta. Using (C5), we transform (C4) as follows: 

AU ^ AW" = --^ imE -mB + P^) - + + -^r, C6 

where E is the bound state energy of quarks which can be obtained from the Schrodinger 
equation with the Hamiltonian Hq. In order to derive (C6) we changed the operator 
by its nonrelativistic expression: p^\E'o = m[E + ^ — Ar — B]^q. As a result of such 
transformations the potential AU2 takes the following form in the case of S-states: 

J/\E'-meson: 

= ^4^-''-' - ^(,nE - ,nB + |^ + + (C7) 

?7c-meson: 



AU2[r) = - e - -—-[mE - mB + 15 ) - -+ 2 ^ + ^ • ^^^i 

A similar transformation of the Breit Hamiltonian can be done for the P-wave states. In 
Table III we present the results of the calculation of the charmonium mass spectrum and 
a comparison with the existing experimental data. The obtained masses agree with the 
experimental ones within an accuracy 1 -h 2 per cent. So, we can suppose that the obtained 
effective Hamiltonian allows to account relativistic corrections in the bound state wave 
functions with sufficiently good accuracy. 
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